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Abstract 
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ON , In this article the particle creation process of scalar and spin 1/2 particles 

in a spatially open cosmological model associated with a universe filled by 
radiation and dustlike matter. The Klein-Gordon and Dirac equations are 
solved via separation of variables. After comparing the in and out vacua, we 

O ' 

Q" 1 obtain that the number of created particles corresponds to Planckian and Fermi- 

Dirac distributions for the scalar and Dirac cases respectively. 
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1 Introduction 



During the last two decades there has been an increasing interest in the study of 
quantum effects in cosmological backgrounds. Since the appearance of the pioneer 
works by Parker and Grib et al among others, a great body of papers have 
been published on the problem of scalar and spin 1/2 particle creation in different 
cosmological configurations. More recently there have been some attempts to analyze 
in this framework the problem of particle creation in the vicinity of cosmic strings 
II 

Since the problem of particle creation in curved space-times is closely related to 
the definition of particle in a curved background, different approaches and interpre- 
tations have appeared during these last two decades, being the most popular the 
adiabatic and the corpuscle interpretations, an overview of the most important re- 
sults and approaches can be found in || and @. Despite the great effort devoted to 
a better understanding of quantum phenomena in expanding cosmological universes, 
the number of articles discussing concrete models is relatively scarce. This is due to 
the fact that, in most of the cases, the relativistic wave equations in curved back- 
grounds are very difficult to solve and therefore is necessary to reduce the original 
problem to a simpler one. 

Recently, the exact solution of Klein-Gordon equation in an expanding (asymp- 
totically flat) cosmological universe filled with radiation has been reported |7J|, the 
authors obtain a thermal planckian distribution of particles created when they com- 
pare the "in" and "out" vacua via the Bogoliubov coefficients. The model presented 
in J7[ encourages us to analyze more complex configurations. 

A good scenario for discussing the problem of scalar and spin 1/2 particle creation 
in a cosmological background, is an isotropic universe with zero cosmological constant 
(A = 0) filled with a mixture of radiation and dust-like matter M, [J in a spatially 
open Friedmann Robertson- Walker model of the form 

ds 2 = -dt 2 + a 2 (t) (d X 2 + sinh 2 (x)(*? 2 + sin 2 ^ 2 )) (1-1) 
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from (|1 . 1|) we obtain 



de 
p + e 



K€ 



-3- 



, da. 



da 
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(1.3) 



where k is the constant appearing in the Einstein's Equations, e is the energy, and p 
is the pressure. Assuming that there is not interaction between radiation and dust, 
we have e = + e r where, H 



A 



B 



a° a H 

after introducing the conformal time r\ related to t by 

dt 



;i.4) 
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and integrating a(r]), and imposing as initial condition a(f]) — > 6 when 77 
being 6 a constant, we obtain the line element 



;i.5) 



-oo, 



c/s 2 = b 2 ^ + l) 2 ( — c/r7 2 + d X 2 +sinh 2 x(^ 2 + sin 2 ^ dip 2 )) (1.6) 

Among the advantages of the line metric ( |1.6| ) we can mention that, this line element is 
regular for any value of the time parameter r]. . The homogeneous and isotropic space- 
time ( |1.6|) evolves in permanent expansion towards and asymptotically flat region at 
rj — > +oo. This asymptotic behavior resembles the one obtained in the perfect fluid 
cosmological model 0. In addition, d/drj approaches a timelike Killing vector when 
f] — > — oo. Also, we have that, after separating variables in the Klein Gordon and 
Dirac equations, the equations governing the time dependence are solvable in terms 
of special functions. It is the purpose of the present paper to obtain the quantum 
spectral distribution of scalar and spin 1/2 particles created in the cosmological model 
associated with the line element (|1.6|) The paper is organized as follows: In §2 we 



solve the Klein-Gordon equation in the background field ( |1.6|) and, in §3 we also solve 



the Dirac equation in this metric. Finally, in §4, based on the results obtained in §2 
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and §3, we compute, via the Bogoliubov coefficients, the probability distribution of 
the number of created particles 



2 Solution of the Klein-Gordon equation 

The covariant generalization of the Klein-Gordon equation reads: 

( 9 aP V a Vp - £R - m 2 ) $ = (2.1) 
where R is the scalar curvature 

R = V(ev\ 1)3 (2 ' 2) 
and £ is a dimensionless parameter which for a conformally coupling takes the value 
off = 1/6. 

Introducing a solution of the form 



$ = T{ V )X(x)Y(&,<p) (2.3) 
and substituting Q2.3| ) into, ( |2.1| ) we find that, after separating variables, Eq. (f2 



reduces to the following set of ordinary differential equations 



d 2 d 1 d 2 \ 



d'd 2 dfi sin 2, & dip 2 

(<P_ + 2 ad + a + ^ + (j2 _ 1 \ T = (25) 
V dr] z adr] a J 

( ' p +2cothx jL + l l±R + A x = (2 . 6) 



\dx 2 dx sinh 2 x 

where I and a are constants of separation, with I > and o 2 > 1. 
After introducing the auxiliary function Z(rj), 



T M = M (2 . 7) 
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^ + fflV(i)) + (i 2 -l\Z = {21 



Eq. ( p.5| ) takes the form: 
substituting the expansion factor a{rf) = b(e v + 1) into (|2.8| ), we obtain the equation, 

m 2 6 2 e 2, + 2m 2 b 2 eV + b 2 m 2 + 2 _ \ % = q 

\d7] 2 J 
whose solution can be expressed in terms of confluent hypergeometric functions 
M(a, b, z) and U(a, b, z) |TT[ as follows: 



Z(r]) = e' ?c exp(-im6e r? )(ciM(^ + c + ibm, 2c + 1, 2ibme n ) 

+ c 2 U{^ + c +ibm,2c+l, 2ibme v ) (2.9) 
where c\ and c 2 are arbitrary constants and c is given by the expression 



c = -z( a 2 + mV - 1) 1/2 (2.10) 

The solutions of (2.4) are given in terms of the harmonic polynomials: Y = Y/ m ($, ip), 
also we have that the solution of the equation (2.6) reads 



x = 7^ P ^ imshx) (211) 



where P^(coshx) are the torus functions [12]. and v 2 = a 2 -I. 



3 Solution of the Dirac equation 

The covariant generalization of the Dirac equation in curved space-time reads |T^| 

{fX^-rg+m}^ = (3.1) 
where 7 M are the curved gamma matrices satisfying the relation, 
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{f\ri + = 2<r (3.2) 

and are the spin connections | I3fl . Choosing to work in the diagonal tetrad gauge 
for 7 M , 

7°= a-V, f = o-V,7 2 =a- 1 (sinh X )- 1 7 2 , (3.3) 
7 3 = a _1 (sinhx) _1 (sin^) _1 7 3 (3.4) 

where 7 M are the standard Dirac flat matrices, we have that the Dirac equation (3.1) 
takes the form 



i - 7 °9„ + -7^ x + — 1— 7 2 ^ + — -1— - 7 3 ^ + m W = (3.5) 
[a a a smh x a smh x sin i> J 

where we have introduced the spinor 



* = a" 2 (sinhx) _1 (sin^)" 1/2 ^ (3.6) 
Applying the algebraic method of separation of variables [0, |I3|, Tfi\ it is possible to 



write eq. ( |3.5D as a sum of two first order differential operators Ki K2 satisfying the 
relation 

[K 1 ,K 2 ]_ = 0, {K! + K 2 }$ = (3.7) 

Ki$ = A$ = -K 2 $ (3.8) 

where 

$ = 7V7 3 $, (3.9) 

K x = (7% + am) 7W, (3.10) 

K 2 = ( 7 ^ x + — ^-7 2 <9, + . 1 . 9 7 3 ^)7S 2 7 3 (3.H) 
sinnx smh x sin w 
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in this way, we have separated the time variable rj from the spatial $j an d f 
variables. The problem arises when we try to reduce ( |3.11| ) to the form ( |3.8j ). It is 
not difficult to see that there are not separating matrices allowing that step. In order 
to go further we write equation (K2 + A)$ = as follows: 



( K 3 + K 4 7 1 7 2 ) $ = (3.12) 
where K3 and K4 are two commuting differential operators given by the expressions 



K 3 = (sinh x) ( 7 2 7 3 <9 X + A ) , K 4 = ( 7 2 7 3 ^ + -?—d v , . 
After introducing the auxiliary spinor S 



(3.13) 



$ = (K 4 + 7 VK 3 )S 
we reduce equation ( |3.12|) to the form; 



(3.14) 



with 



where 



Ni + N 2 E = 0, Ni, N 



N 2 E = -NiE = k 2 E 



(3.15) 



(3.16) 



Ni = (sinhx) 2 ( 5 2 + cothx<9 x — A 7 2 7 3 cothx + A 2 ) 

N 2 = (9| + z7 2 7 3 fc (/ ,(sin^)- 2 cos?9-fc 2 (sin^)" 2 ) 

ktp in Q3.18 ) is the eigenvalue of the operator — id^. 

Choosing to work in the following representation of the Dirac matrices: 

/ 



7 



V 
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(3.17) 
(3.18) 
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we have that the spinor S takes the form 

s = q(x)B(#) 

the functions p(x), <l{x)i r (x) an d satisfy the differential equation 



(3.19) 



(3.20) 



(smhx(d x + Ti\))(smhx(d x + ±t\))X(x) + k 2 X(x) = (3.21) 

where the upper signs correspond to p(x)) r (x); an d the lower signs are related to 
q(x) and s(x). 

The functions B ($), C($), are solution of the equation 



(d 1? ±^(sin^)- 1 )(^^^(sin^)- 1 )r(^)-A; 2 r(^) = (3.22) 

where the upper signs correspond to A(d) and C($), and the lower ones correspond 
to B($) and D(-&). After substituting the expression Q3.20|) into (|3.14j ) we arrive at 



ik 



(3.23) 



q( X )A('&)(l-i)a 

where a and (3 are functions depending on the variable rj. The functions A, B, C, D, 
and p, q, r, s and are related by the following systems of coupled partial equations: 



sini?' 



C 




k ( B 

(do + At) 

smv \ D 



C 



(3.24) 
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sinh X {d x + iX) \ P ) = fc ( J , sinh x(d x - iX) ( ] = -jfe f P j (3.25) 
r / \ s / \ s I \ r 



where e and to are constants. From ( |3.23|) and (|3.10|) we can put p(x) = 7*(x) 



<z(x) = ^4(^) = C(^), and £>(■#) = D($). The functions a and /3 satisfy the 
system: 

(d. n + ima)f3 = — Xa, (d n — ima)a = X [3 (3.26) 
from which we obtain: 

[d\ + m W 7 + e" (2mV - z'mfc) + m 2 b 2 + A 2 )a = (3.27) 

(d 2 v + m 2 b 2 e 2r > + e v (2m 2 b 2 + imb) + m 2 b 2 + X 2 )f3 = (3.28) 

whose solution can be written in terms of confluent hypergeometric functions as fol- 
lows 

a( v ) = e^e-*"**" {-co ^ - ' ^ Af (c + zmfc + 1, 2c + 1, 2im6e'») 

A 

(V 2 _|_ m 2 h 2 ) 

+ ct± J -U{c + imb + 1,2c + 1, 2im6e ,? )) (3.29) 

A 

/5(r/) =e" c e- im6e " (c M(c + im&, 2c + 1, 2imbe v ) + c x U{c + imb, 2c + l,2imbe T >)) 

(3.30) 

where c is given by the expression: 

c = _ 2 (A 2 + m^ 2 ) 1 ^ (3.31) 



The solution of the coupled system of equations ( |3.24 ) can be expressed in terms of 
Jacobi polynomials as follows: 

A{&) = sinV(^)cos(?9/2)P,(V 1 / 2 'V _1/2) (cos^) (3.32) 
B(&) = sinV(^)sin(^/2)PiV- 1/2 'V +1/2) (cos^) (3.33) 
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where n is given by the expression 



n 



The system (|3~25p can be rewritten in the following matrix form 



(3.34) 



I oL + iXcr" 



i k 



V 



sinh x 

after applying the matrix transformation S 



-a' 



(3.35) 



s- 1 



V 



p 



1 



we reduce ( |3.25| ) to the form 



ik 



sinh x 



)Q 



S = —(l+ia 1 
Ql V2 



XP (d x + ~^-)P 
v x smhx 



XQ 



with, 



(3.36) 



(3.37) 



P + iQ= p, Q +iP 



(3.38) 



the solution of the system (|3.37| ) can be expressed in terms of Gauss hypergeometric 
functions [ P2 |, 



P( X ) = sinh lk ( X ) cosh ( x /2) x 

113 

F(- + iX +ik,- -iX + ik,- + ik, (I -coshx)/2) 



Q(x) 



a + ik)/\ ) sinh ik { X ) sinh(x/2) x 



+ iX +ik,^ -iX +ik,^ + ik, (I -coshx)/2) 



(3.39) 



(3.40) 



Finally, from ( |3.6j ), ( |3.9| ) and 23| ) we find that the solution of the Dirac equation 
Ah reads 
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= a 2 (sinhx) (sini?) 



-1/2 



(3.41) 



I p(x)B^)(i + l)^( V ) ^ 
q( X )A(#)(l-i)[3{r]) 
p( X )B(-&)(l + i)a( V ) 
\ q(x)A(#)(l - i)a( V ) J 
Here some comments are in order. The election of a diagonal tetrad for separating 



variables and solving the Dirac equation in the background field ( |1.6|) was done for the 
sake of simplicity. If we want to reobtain an angular dependence of the spinor ( |3.41| ) 
in terms of spherical harmonics, it is enough to achieve the transformation relating 
the diagonal tetrad gauge with the Cartesian one |17| Obviously, this transformation 
does not affect the time dependence of the spinor solution. 



4 Particle creation 

In this section we proceed to analyze the phenomena of scalar and Dirac particle 
creation in the cosmological background associated with the metric ( |1.6|) . The model 
in question begins at rj = — oo and evolves to rj — > +oo. Among the particularities of 
our model we have to mention that, like @, the universe (|i~6l ) is asymptotically flat, 
but does not present a period of contraction during its evolution toward the future. 

In order to establish the asymptotic behavior of the solutions, at rj = — oo and 
at 7] = +oo, of the corresponding wave equations, we write the Hamilton- Jacobi 
equation 

g^S, a S, p + m 2 = (4.1) 



in the background field (1.6). Then we obtain: 



-(S, v ) 2 + (S, x ) 2 + — — US,,) 2 + — (S^) 2 )+ m 2 a 2 = (4.2) 
We can separate in Eq. ( |4.2|) the time dependence from the spatial one as follows: 
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then, we obtain 



S = Tfo) + F{x,#,<p) 



(4.3) 



(T,,) 2 = mV + A 2 , — > T = ± /(mV + A 2 ) 1/2 dr] (4.4) 

where A is a constant of separation. Then, the quasiclassical asymptotes of the 
solutions are 



* -> F(x, 0, exp ( ±z /(m 2 a 2 + A 2 ) 1/2 dr] ) (4.5) 
The asymptotic behavior of the function T(?7) is 

lim T -> ± (m 2 6 2 + A 2 ) 1/2 r?, lim T -> im&e'' (4.6) 

7/— > — OO IJ— ! >00 

The relation between the conformal time 77 and the cosmic time t can be obtained 
from ( |1.5|) giving as result, 

b{e* + v ) = t (4.7) 

notice that the functions b(e v + rj) monotonically increases, and t — > ±oo when 
7] — > ±oo, and, when 77 — > +oo the metric ( |1.6|) becomes a Milne universe || 

We have that the solutions of the Klein-Gordon equation with the asymptotic 
behavior (|4.5| )can be written as, 



Z^rj) = (2c)~ 1/ V c exp(-im&e , ')M(i + c - ibm, 2c + 1, 2ibme n ) (4. 
for Tj — > — oo, and 



Zoo(v) = e^ c/ V c exp(-im&e r? )?7(- + c - ibm, 2c + 1, 2i6me ?? ) (4.9) 



for fj — ► oo. where, in order to analyze the form ( |4.8|) takes at infinity, we have used 



the formula: [ I I 
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U(a, 7, z)^z- a (l + 0{ l -)y (4.10) 

and the relation (|4.7|) 

Notice that, for the constant c given by (p . 1 Of) , T_oo and can be associated with a 
positive frequency mode at i] — > — oo and 77 — > oo respectively. The next step is to 
compute the Bogoliubov coefficients relating the two different vacua at 77 = — 00 and 
i] — > 00 m . In the present case it is possible to obtain (3 without performing the inte- 
gration. Indeed, using the recurrence formula between the confluent hypergeometric 
functions U(a, 7, z) and M(a, 7, z) [i ll]. 



M(a, 7 ,z) = r(7) e^Ufa^z) + ^M^-aj,-?) (4.11) 
1 (7 — ctj 1 (aj 



substituting a = \ + c +ibm, 7 = 2c + 1 and 2 = 2ibme v into ( 4.10|) and 
taking into account the relation 

e z U(<y -a,j,-z) = eV -7 (-l) 1_7 l7(l - a, 2 - 7, -z) (4.12) 

we arrive at, 

+ (2c) v4 2 + + ^ -^^MXW (4-13) 



then, using the relations for the T function, [12 



2 2z - 1 T(z)T(z + \) = V^T(2z) (4.14) 



T (W) \ 2 = ■ ^ S (4-15) 



and considering the limit case when b is small, we arrive at 
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i p i 2 = ( 4 - 16 ) 

where is u — (a 2 — l) 1 / 2 

Analogously, we have that the Bogoliubov coefficients for the Dirac case can be ob- 
tained by means of the functions a(rf) and [3 (if) determining the time dependence of 
the spinor (|3.41|). Then, considering the corresponding positive frequency solutions 



at the origin and at infinity, 

£-00(77) = e r ' c e- imbeV M(c + imb,2c+l,2imbe r ') (4.17) 
Pociv) = e^e^e-^Uic + imb^c+l^imbe 71 ) (4.18) 
where the constant c is given by ( |3.31| ). Taking into account the relation Q4.10 ) and 



( |1.11| ) we obtain 



Pa(v) = J^e'^e-^M + 
1 (c + 1 + imb) 

+ r ( 2c + U_ e-^ c + mb +^e inc/2 (2mb)- 2c f3'(r]) (4.19) 
T(c — imb) 

then, for small values of b, using the relation for the gamma function T(z), |12| 

|r(^ + i)| 2 =— ^- (4.20) 

2 cosh(7n/J 



and the expression ft4.14|) we arrive at: 



P I 2 = (4-21) 



1 

2 e 27rA + 1 

where A in ( f4.21| ) is the separation constant appearing in ( |3~8| ) 
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5 Concluding remarks 



In this paper we have analyzed the process of creation of scalar and Dirac particles in 
an isotropic and homogeneous cosmological universe filled by non interacting radiation 
and dust-like matter. It is worth noticing that, also in this model we have obtained 
that the average number of particles created satisfies a thermal distribution law. This 
result show that the equation of state does not determine the quantum spectrum 
but the time behavior of the line element associated with the cosmological universe. 
From a mathematical point of view, this paper shows the capabilities of the algebraic 
method of separation of variables and how it works when first order separation is 
not possible. We hope that the results presented in this article will be of help in 
understanding quantum phenomena in cosmological backgrounds. 

References 

[1] Parker L. 1968 Phys. Rev. Lett. 21, 562 

[2] Grib A. A., Mamaev G. G., and Mostepanenko. V. M. (1970) Soviet J. Nuclear. 
Phys. 10, 722. 

[3] Parker L. 1989 Phys. Rev. Lett. 59, 1369 

[4] Sahni V. 1988 Mod. Phys. Lett. 3, 1425 

[5] Birrel N. B. and Davies P. C. W. 1982 Quantum Fields in Curved Space (Cam- 
bridge: Cambridge University Press) 

[6] Grib A. A., Mamaev G. G., and Mostepanenko. V. M. Quantum Vacuum effects 
in strong fields (Energoatomizdat, Moscow, 1988) 

[7] Costa I. Deruelle N. Novello M. Svaiter N. 1989 Class. Quantum. Grav. 6 1893 



14 



[8] Chernin A. D. 1966 Sov. Astron. 9, 871 
[9] Kharbediya L. I. 1977 Sov. Astron. 21, 672 
[10] Parker L. 1975 Phys. Rev. D. 12 1519 

[11] Abramowitz M. and Stegun I. A. Handbook of Mathematical functions (Dover, 
New York, 1964) 

[12] Gradshtein I. S. Ryzhik I. M. Table of integrals, Series and Products (Academic 
Press, New York 1980) 

[13] Brill D and Wheeler J A 1957 Rev. Mod. Phys. 29 465. 

[14] Andrushkevich I. E and Shishkin G. V 1987 Theor. Math. Phys. 70 204 

[15] Shishkin G V and Villalba V M 1989 J. Math. Phys. 30 2132 

[16] Shishkin G V and Villalba V M 1992 J. Math. Phys. 33 2093 

[17] Villalba V. M. and Percoco U. 1990 J. Math. Phys. 31 715 



15 



